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Abstract 

Logan's graphical analysis (LGA) is a widely-used ap- 
proach for quantification of biochemical and physiological 
processes from Positron emission tomography (PET) im- 
age data. A well-noted problem associated with the LGA 
method is the bias in the estimated parameters. We re- 
cently systematically evaluated the bias associated with the 
linear model approximation and developed an alternative 
to minimize the bias due to model error. In this study, 
we examined the noise structure in the equations defining 
linear quantification methods, including LGA. The noise 
structure conflicts with the conditions given by the Gauss- 
Markov theorem for the least squares (LS) solution to gen- 
erate the best linear unbiased estimator. By carefully taking 
care of the data error structure, we propose to use struc- 
tured total least squares (STLS) to obtain the solution us- 
ing a one-dimensional optimization problem. Simulations 
of PET data for [11C] benzothiazole -aniline (Pittsburgh 
Compound-B [PIB]) show that the proposed method signif- 
icantly reduces the bias. We conclude that the bias associ- 
ated with noise is primarily due to the unusual structure of 
he correlated noise and it can be reduced with the proposed 
STLS method. 



1. Introduction 

Graphical analysis (GA) is a routine tool for quantita- 
tive imaging with PET in various clinical and physiological 
studies. The first GA method is the commonly usedd Patlak 
method that was introduced by Patlak, ifTTl [lOl . for irre- 
versible tracers. Logan extended this method for reversible 
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tracers, Q. So far a set of GA methods have been devel- 
oped for both reversible and irreversible systems, for plasma 
input and reference model, and for calculation of uptake 
rate, distribution of volume (DV) and DV ratio (DVR) or 
binding potential (BP). 



A well-noted problem with the use of GA methods, par- 
ticularly for reversible system, is the bias in the estimated 
parameters, 0131 151, To reduce the bias, Logan et al suggest 
to deal with noise by smoothing the data, J6), and Varga et 
al proposed perpendicular least squares, 03), which is ex- 
actly the total least squares (TLS) method, (2). Ichise et 
al, J5), rearranged the equation to a multilinear equation to 
decrease the noise, (5). Ogden used a nonlinear likelihood 
estimation, 0. The bias associated with GA approaches, 
we believe, has three possible sources. Published work pri- 
marily dealt with the bias related to the random noise, one 
of the three sources. The other two sources are the nu- 
merical quadrature error and an approximation of the un- 
derlying compartmental model. We recently systematically 
evaluated the bias associated with the model approximation 
and developed an alternative model for minimizing the bias 
caused by model error, ||3|. In this study, we investigate the 
noise effects in parameter estimation differently. In con- 
trast to the linear least squares (LLS) or ordinary TLS algo- 
rithms, careful examination of the data error structure leads 
to our proposal to develop a structured total least squares 
(STLS) approach to estimate the parameters. Simulation 
shows that the bias due to noise is greatly reduced by the 
STLS method. 



The rest of the paper is organized as follows: The new 
approach, STLS, is introduced in Section |2j The simula- 
tion study is described in Section [3] and results reported in 
Section |4] Issues relevant to the proposed approach are dis- 
cussed in Section|5] Conclusions are presented in Section|6] 
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where I Pi = J Q * C p (s)ds. In matrix notation this is 



2.1 Logan's method and alternative linear 
methods 

Assume the equilibrium is reached at some time point 
t', after which linear equations associated with corre- 
sponding linear methods are assumed valid. Logan's GA 
(LGA) quantification method for reversible radiotracers 
with known plasma input function, i.e. plasma concentra- 
tion of the unmetabolized tracer, is based on the following 
equation 



LGA 



L Ct(s)cIs L C v (s)ds 

~ V „ ,. , h b, (1) 



Or(t) C T (t) 

where Cr(t) is the measured "tissue time activity curve" 
(TTAC), C p (t) is the input function, and V represents for 
the volume distribution (DV). The equation approximately 
reflects the tracer behaviors over the equilibrium period. 
Ichise et al revised (Q]i to a multilinear equation as follows , 
0, 

MAI :C T (t)^-jj o C p (r)dr + i J C t {t)At. (2) 

Both LGA and MAI are solved by LLS. In the ordinary 
least squares method the independent variables are assumed 
to be noise-free. But Varga et al, fl4l . noted that noise ap- 
pears in both the independent and dependent variables and 
proposed to use TLS for equation (01, let us call it MAO, 
which is also the root equation for LGA and MAI: 

MAO: f CT(T)dr«Dv/ C p {r)dT + bC T (t). (3) 
Jo Jo 

The TLS solution is obtained from the right singular vector 

corresponding to the smallest singular value of the matrix 

formed by columns f Q Ct{t)&t, J* C p (t)(1t and Ct(£). 

A complete introduction and analysis of basic algorithms 

for TLS is presented in (4) . 

2.2 Method development 



We denote the durations and central time of the scanning 
frames by Aj and ti,i — 1, • • ■ ,n and assume t' = t q . 
Thus ti,i = q, ■ ■ ■ ,n fall in the equilibrium period. Be- 
cause Ct(U) = Ct(s)^s / Ai the integral 7C, = 

J U C T (s)ds = EjpiOrfeOA,- + C T (U) Ai/2. With 
these notations the discretized multilinear equation of (0 
can be written as: 



f IP q , C T (t q ) ^ 
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To simplify the discussion we assume 1) the equilibrium 
is truly reached after t', thus the discussed linear models 
do not have model error; and 2) the noise in C p (t) can be 
ignored as compared with the noise in Cr(i). The Gauss- 
Markov theorem tells us that the LLS estimator is the best 
linear unbiased estimator assuming that no noise in the in- 
dependent variables and the entries of the noise vector are 
normally distributed i.i.d. variables with zero mean and 
common variance. Because both Cr(i) and JjJ Cr(s)ds 
contain noise and the noise is correlated, none of Gauss- 
Markov theorem's conditions are satisfied for the data in 
linear equations (dJ, (f2]i and ([3]). This explains the bias of 
these linear methods. On the other hand, the TLS solution 
is equivalent to the maximum likelihood solution when the 
noise in independent and dependent variables are i.i.d mul- 
tivariate normally distributed with zero mean and common 
covariance, [;4|. The noise in the linear equations do not 
satisfy this condition either. Thus TLS does not produce an 
unbiased parameter either, as noticed in 0. 

Let Ct(U) + fi be the true radioactivity at time ij, i.e. 
—fi is the measurement error in Ct(U)> assumed to be ap- 
proximately normal distribution with zero mean and vari- 
ance var(fi) = af. Thus 7 = J^i 1 fi^i nas variance 
X^i _1 ( f7 «Ai) 2 . By incorporating the noise in the approxi- 
mated equation (0]i we obtain the following exact equation, 
which reflects the structure of the noise, 



A[V,b] T + M = d + L{ + 7 1, 



(5) 



where f = [f q , ■ • • , f n ] T , 1 = [1, ■ • ■ , 1] T , and L is a lower 
triangular matrix: 

/ A,/2, 



L = 



A, 



A,+i/2 
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We propose to estimate V by solving an optimization prob- 
lem as follows: 



STLS 



mm w"f 
subject to 



-{ T Ei 



(6) 



where E = Aiag(\ f a 2 q , ■ ■ ■ and w = 

l/(Si~ 1 ( cr iAi) 2 ). © is a STLS problem and equivalent 
to maximum likelihood (ML) model, as demonstrated 
in section [5] The n — q + 4 unknowns are V, b, 7 and 



fi,i = q,- ■ ■ ,n. For example, if 3 frames fall in the equi- 
librium period, i.e. n — q + 1 = 3, we have 6 unknowns. If 
we fix the variable b the optimization problem (|6]l becomes 
a quadratic programming problem and can be simplified to 
a LLS problem as follows. Let us use ai and a 2 to denote 
the two columns of matrix A, i.e. A = [ai,a2]. Equation 
d5} can be rewritten as 

(L-bl)i =[a u -l][V, 1 ] T + a 2 b-d. 



LetB = (L-bI)- l [aii,-l] andg 
the objective function becomes 



(L-bI)- 1 (a 2 b-d) 



F(b,V,j) 
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By setting the first order derivatives of F(b, V, 7) with re- 
spect to V and 7 to zero we obtain the following equation 



(B T EB + wJ) 



where matrix 



J 



V 




1 



-B T Eg, 



(7) 



Denoting the solution of © by V(b) and 7(6) we define 
function G(b) = F(b,V(b),j(b)). Based on the above 
analysis we design a numerical algorithm for solving the 
STLS problem (O by an one dimensional minimization as 
follows 

Algorithm 1 Given b e [Pi, fa]. 



1. Solve 



min G{b), for b G [fa, fa]. 



2. For the solution of the above minimization, b*, calcu- 
late corresponding V by (j7J. 



only exception is the ROI 6, which is far from equilibrium 
even after 200 minutes. For the noise-free decay-corrected 
concentration TTAC, CJ(i), Gaussian noise at each time 
point tj, G(0, cr(C^{t)), is modeled using the approach in 
@[I1|5] as follows 



a(Ci(ti)) = Sc x 



IC*{U)e 
60AU 



(8) 



here A is the tracer decay constant (0.034 for 11 C) and Sc 
is a scale factor, which is set to 0.5 and 1 in our simulations. 
The units for Cx(ij) and ti are kBq/ml and minutes respec- 
tively and Atj is multiplicated by 60 to measure the frame 
duration in seconds. 1000 random sample sets are tested for 
each noise level, Sc=0.5 and Sc=l. 

4 Results 

We present histograms for the percentage relative error 
of the bias, 100(V est - W) /DV, in FigureQ]with the range 
of the percentage error for each method indicated in the leg- 
end. The upper figure is for noise scale Sc = 0.5 while the 
lower figure is for Sc = 1. It is clear that the variances for 
the results of MAI and TLS and the bias of the results of 
LGA are too large. STLS outperforms all these three meth- 
ods. There are some situations, in which the relative error 
is less than —100%; in other words, the calculated DVs are 
negative. This occurs 133 and 29 times over all 11000 ran- 
dom tests for MAI and TLS at noise scale Sc = 0.5 while 
509 and 95 occurrences are observed for MAI and TLS for 
noise scale Sc = 1. LGA and STLS do not produce any 
negative DV. The reason for the negative DV for MAI is 
discussed in (3J- In the simulations, the average CPU time, 
in seconds, per TTAC were 4e — 4, 2e — 4, 3e — 4 and 
9.4e - 3, for LGA, MAI, TLS and STLS, respectively. 



3 Simulation study 



5 Discussion 



The simulated data are adopted from published clinical 
data, I1T21 [TBI . Specifically, arterial input function and rate 
constants of two tissue reversible compartmental model for 
PIB tracer are used. Eleven regions , ROI 1, to ROI 11, of 
normal controls (NC) and Alzheimer's Disease (AD) diag- 
nosed subjects are tested. Details of the data are described 
in@. 

Given the decay corrected input function and the kinetic 
parameters we generate corresponding unperturbed TTACs, 
C^(i). The frame durations, total 240 minutes scanning, are 
set to be, given in minutes, 4 x 0.25, 8 x 0.5, 9x1, 2x3, 8x5 
and 18 x 10. Frames falling in 120 to 240 minutes, i.e. frame 
37 to 49, are chosen as equilibrium frames. We set t' = 
120 so that the equilibrium is approximately attained. The 



The proposed STLS is equivalent to the ML solution 
with the Gaussian density assumption for noise /j. Because 
the variances for /,; and 7 are of and E 2 = J^i 1 (c r iAj) 2 
respectively, the ML problem can be formated as follows 

max e T 2 /£ 2 e Er=,A 2 /^ 2 
subject to (J5). 

The equivalence is easily shown by taking the logarithm of 
the objective function. 

In our simulation the exact variances are used. In prac- 
tice, the variance of each frame can be estimated, H] |9) . In 
this work we assume that the linear simplification for MAO 
does not have significant model error, which is the reseason 
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Figure 1. Histograms for normalized error (in 
percentage),i00(V cst -DV)/DV, of the results 
for all eleven ROIs and four methods, LGA, 
MA1 , TLS and STLS. The upper and lower fig- 
ures are corresponding to noise scale Sc = 
0.5 and l respectively. The error ranges are 
presented in the legends. The errors are as- 
signed to 7 bins with centers specified at 
-80%, -50%, -20%, 0%, 20%, 50% and 80%. 
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Figure 2. Function G(b) for a representative 
case. 



we use a long scan duration 240 minutes. If the model error 
does exist, this is the case for practical situations, we need to 
correct the model error, (3), and perform further research to 
reduce the noise effects based on the model error corrected 
model. At last, in our simulations the range for b is given by 
[—150, —10] for all cases. A representative curve for func- 
tion G(b) is illustrated in Figure[2] which is a convex curve. 
If a case dependent good range of b can be estimated the 
performace of the STLS algorithm can be further improved. 

6 Conclusions 

In this article we proposed to reduce the noise effects of 
linear PET quantification by a structured total least squares 
noise model and developed an efficient numerical algorithm 
for its solution. We validated our findings through simula- 
tions with clinical derived PIB-PET data. Simulation re- 
sults demonstrate that the STLS algorithm significantly re- 
duces the bias caused by noise in PET data as compared 
with LGA, MAI and TLS. 
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